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Abstract 

. Quantum-mechanical multiple-well oscillators exhibit curious com- 

plex eigenvalues that resemble resonances in models with continuum 

•i— I . 

^ ' spectra. We discuss a method for the accurate calculation of their real 

H ; 

5t , s-iid imaginary parts. 



1 Introduction 

Some time ago, Benassi et al [1] discussed the occurrence of complex eigen- 
values, or "resonances ", in some quantum-mechanical multiple-well oscilla- 
tors, and calculated them for a particular example. Recently Killingbeck [2] 
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showed that the Hill-series method yields quite accurate results for both the 
real and imaginary parts of those eigenvalues if one introduces a complex pa- 
rameter in the exponential factor of the expansion. In principle, one has to 
tune up this parameter in order to obtain an acceptable rate of convergence. 
Such " complexification" of the well-known Hill-series method had been tried 
successfully before in perturbation and matrix approaches [3-5]. Complex- 
ification is a term coined to indicate the use of, for example, a complex 
frequency in the treatment of a perturbed harmonic oscillator or a complex 
atomic number in the case of a perturbed Coulomb problem [2-5]. 

The Riccati-Pade method (RPM) is known to be suitable for the accurate 
calculation of bound states and resonances of simple quantum-mechanical 
models [6-14]. However, it has only been applied to the most commonplace 
resonances in the continuum spectrum [10-14]. The purpose of this paper is 
to investigate if the RPM is also a reasonable alternative to the calculation of 
the unusual kind of resonances considered by Benassi et al [1] and Killinbeck 
[2]. 

In Section [2] we outline the RPM and in Section [3] we apply it to the three- 
well oscillator treated explicitly by Benassi et al [1] and Killingbeck [2]. 

2 The Riccati-Pade method (RPM) 

In order to make this paper reasonably self-contained, in this section we 
outline the RPM in a quite general way. Suppose that a solution to the 
eigenvalue equation 

^"{x) + [E- V{x)] ^{x) = (1) 
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can be expanded in the form 

oo 

^(a;) = ^ CjX^^ a,p>0. (2) 
The power-series expansion for the regularized logarithmic derivative 

converges in a neighbourhood of x = and the coefficients fj depend on the 
eigenvalue E. The function /(x) is a solution to the Riccati equation 

/'(x) - f{xf + -/(x) + V^(x) -E~ <^i^_zA = 0. (4) 

Equations ([I])-(j4]) apply to both one-dimensional (— oo < x < oo) and 
central-field (0 < x < oo) models. If V{x) is a parity-invariant one- 
dimensional potential, then a = for even states, a = 1 for odd ones, 
and /5 = 2 for both cases. If lima.^o+ x'^V{x) = V^2 > 0, then a{a — 1) = V-2 
removes the singularity at origin in the case of a central-field model. 

The RPM consists of rewriting the partial sums of the power series ([3]) 
as Fade approximants x^~^[N + d/N]{z), z = x^, in such a way that 

-r^N+d j 2N+d+l 

[N + d/N]{z) = ^;\^ ^ = E fjz^ + 0{z'^^'^^). (5) 

In order to satisfy this condition the Hankel determinant H'l^, with matrix 
elements fi+j+d+i, hj = 0, 1, . . . , iV, vanishes, where D = + 1 = 2, 3, . . . 
is the determinant dimension, and c? = 0, 1, ... is a displacement [6-14]. The 
main assumption of the RPM is that there is a sequence of roots El^-'^l of 
the Hankel determinants -fT^ that converges towards a given eigenvalue of 
the Schrodinger equation as D increases [6-14]. For brevity we call it a 
Hankel sequence. 
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Notice that one obtains the coefficients fj from the expansion of the 
Schrodinger equation ([1]) or the Riccati equation (jl]) quite easily, and that 
unhke the Hill-series method [2] the RPM does not require an adjustable 
complex parameter. Besides, it is not necessary to take into account the 
boundary conditions explicitly in order to apply the RPM, and, for that 
reason, the method provides both bound states and resonances simultane- 
ously [6-14]. 

3 Results and discussion 

In what follows we apply the RPM to calculate the curious complex eigen- 
value of the triple-well oscillator 

V{x) =x^ - Ig^x^ + g^x^ (6) 

reported by Benassi et al [1] and Killingbeck [2]. In this case {3 = 2 and we 
choose a = for even states as discussed above. 

Table [1] shows a Hankel sequence E^^''^^ that converges towards the low- 
est complex eigenvalue when g = 0.14. We have kept twenty digits in all 
entries in order to show how they become stable as D increases. Notice the 
remarkable rate of convergence of the Hankel sequence for both the real and 
imaginary parts of the eigenvalue. 

Table [2] shows the same complex eigenvalue for a range of (^-values some- 
what wider than the ones chosen by Benassi et al [1] and Killingbeck [2]. We 
have truncated the results, obtained from Hankel determinants with D < 15 
and d = 0, to the apparently last stable digit. The first digits of our results 
agree with those given by Benassi et al [1] and Killingbeck [2]. We notice 
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that lmE{g'^)g'^ exp{l/{2g'^)) does not seem to approach a constant for those 
values of g. It may be that lmE{g^) attains the WKB asymptotics [1] at 
smaller values of g. 

It is interesting to compare the strange resonance of the potential ([6]) 
with the more commonplace one of the potential 

V2{x) = - 2g^x^ (7) 

that was treated earlier by means of the RPM [10]. Table [3] shows the lowest 
resonance for this model for the same values of g considered before. We 
appreciate that the imaginary part of this resonance is considerably greater 
than the previous one and that it seems to approach the WKB asymptotics 
ImE^^^ig^) = [A/{2ng^)]exp{-l/[3g^]) somewhat faster. 

The results of this paper clearly show that the RPM is suitable for the 
calculation of both real and complex eigenvalues of simple Hamiltonian op- 
erators, even in the case of quite small imaginary parts. We believe that this 
approach is a most useful tool in the numerical investigation of a wide variety 
of eigenvalue problems. Its main advantages are: great rate of convergence 
and simple straightforward application that does not require adjustable pa- 
rameters or explicit consideration of boundary conditions. 
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Table 1: Convergence of a Hankel sequence E^^'^^ towards the lowest complex 
eigenvalue of the oscillator ([6]) with g = 0.14. 



D 


ReE 




ImE 




2 


0.96913474062929793208 









3 


0.96912933030952144688 









4 


0.96912932029284635448 









5 


0.96912932006642961226 


3.678: 


L221743857153252 10" 


-10 


6 


0.96912932002647227146 


3.3990326234127550889 10" 


-10 


7 


0.96912932002710973379 


3.3801038698293392418 10" 


-10 


8 


0.96912932002717289039 


3.379J 


B079586780234680 10" 


-10 


9 


0.96912932002717518442 


3.379f 


B093143407212241 10^ 


-10 


10 


0.96912932002717525409 


3.379f 


S095397280767486 10" 


-10 


11 


0.96912932002717525622 


3.379f 


S095479442123313 10" 


-10 


12 


0.96912932002717525629 


3.379f 


S095481219295624 10" 


-10 


13 


0.96912932002717525629 


3.379J 


B095481219029216 10" 


-10 


14 


0.96912932002717525629 


3.379J 


B095481216587093 10" 


-10 


15 


0.96912932002717525629 


3.379J 


B095481216435223 10" 


-10 
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Table 2: Complex eigenvalue of the oscillator ([6]) for several values of g. 



9 




ReE{g ) 




lmE{g ) 


lmE{g^) g'exp{l/{2g^)) 


0.08 


0, 


.99025645954150600314 


1, 


.16994 10-32 


0, 


,6362094894 


0.09 


0, 


.98761765110834730415 


1, 


.28623698 lO-^^ 


0, 


,6700502315 


0.10 


0, 


.98464158830285882643 


1, 


.3513930260 lO^^o 


0, 


,7006574893 


0.12 


0. 


,97763491479323529157 


4, 


.3530125379031 10"^^ 


0, 


,7530467190 


0.14 


0. 


,96912932002717525629 


3, 


.37980954812164 10"^° 


0, 


,7944913345 


0.16 


0, 


,95896997046169207832 


1, 


.0619001732959989 10"^ 


0. 


,8253492417 


0.18 


0. 


,94691604067745932355 


5, 


.18077667159013113 10^*^ 


0, 


,8453084682 


0.20 


0. 


,93255571582477452180 


7, 


.94775543996767651 10"^ 


0, 


,8530716514 


0.22 


0. 


,91525354748034208273 


5, 


.70253065914296141 10"^ 


0. 


,8461088416 


0.24 


0, 


,89442055320991452496 


2, 


.424632840047890532 lO'^ 


0, 


,8222158493 


0.26 


0. 


,87011531157430539225 


7, 


.104058338260953225 10~'^ 


0, 


,7828715436 


0.28 


0. 


,84333442392342060412 


1, 


.5915859465250206010 10^^ 


0, 


,7343132667 


0.30 


0. 


,81560795814733914293 


2, 


.9400216892153485663 10^^ 


0, 


,6844475376 
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Table 3: Lowest resonance of the oscillator ([7]) for several values of g. 



ReE{g'') lmE{g') lmE{g') g exp{l/{3g^)) 



0.08 


0, 


.99017315154568105030 


4, 


.66667951 lO'^^ 




1, 


.554541174 


0.09 


0, 


.98748105548308533216 


2, 


.3014736620 10"^^ 




1, 


.543296673 


0.10 


0, 


,98442766976525540084 


5, 


.1093948883947 10"^^ 




1, 


.530566484 


0.12 


0. 


,97716020191841551216 


1, 


.1063680213861671 10"^ 




1, 


.500354438 


0.14 


0. 


,96816424784205963513 


4, 


.297124100601175228 10"^ 




1, 


.463074727 


0.16 


0. 


,95708500653988706061 


1, 


.9606870293524100682 10" 


5 


1, 


.417112487 


0.18 


0. 


,94328218799381038166 


2, 


.5699864836055797687 10" 


'4 


1, 


.35910675 


0.20 


0. 


,92594246107314318252 


1, 


.5440221243204925966 10" 


3 


1, 


.284707315 


0.22 


0, 


.90482508551985951067 


5, 


.5395017058573660278 10" 


3 


1, 


.193719284 


0.24 


0, 


.88093011197386366807 


1, 


.3978475279423154843 10" 


-2 


1, 


.093828654 


0.26 


0, 


.85613353763295142744 


2, 


.767004146177769213 10"^ 




0, 


.9964939951 


0.28 


0, 


.83225989985769363726 


4, 


.6300611971065823176 10" 


2 


0, 


.9104055713 


0.30 


0, 


.81052712217939364397 


6, 


.8908503646837670242 10" 


-2 


0, 


.839251556 
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